Phase diagram of insulating crystal and quantum Hall states in ABC-stacked trilayer 

graphene 
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In the presence of a perpendicular magnetic field, ABC-stacked trilayer graphene's chiral band 
structure supports a 12-fold degenerate N = Landau level (LL). Along with the valley and spin 
degrees of freedom, the zeroth LL contains additional quantum numbers associated with the LL 
orbital index n = 0, 1,2. Remote inter-layer hopping terms and external potential difference Ag 
between the layers lead to LL splitting by introducing a gap All between the degenerate zero- 
energy triplet LL orbitals. Assuming that the spin and valley degrees of freedom are frozen, we 
study the phase diagram of this system resulting from competition of the single particle LL splitting 
and Coulomb interactions within the Hartree-Fock approximation at integer filling factors. Above 
a critical value All of the external potential difference i,e, for |All| > All, the ground state is 
a uniform quantum Hall state where the electrons occupy the lowest unoccupied LL orbital index. 
For |All| < All (which corresponds to large positive or negative values of As) the uniform QH 
state is unstable to the formation of a crystal state at integer filling factors. This phase transition 
should be characterized by a Hall plateau transition as a function of All at a fixed filling factor. 
We also study the properties of this crystal state and discuss its experimental detection. 

PACS numbers: 73.21.-b,73.22.Gk,72.80.Vp 



I. INTRODUCTION 

Quantum Hall studies of graphene's two dimensional 
electron system provided the earliest confirmation of the 
massless Dirac character of its bandsii^. Similar studies 
have also confirmed the massive Dirac character of bi- 
layer graphene bands'^''' and are making progress toward 
revealing the distinct bands of ABA and ABC trilayers^. 
Few layer graphene systems stacked in different ways 
give rise to distinct electronic properties at relevant en- 
ergy scales which introduces an entirely new class of two- 
dimensional electron gas (2DEG) systems. Recent the- 
oretical work has indicated that except for very weak 
fields, multilayer graphene systemsSi^ can be classified 
as a new class of 2DEG systems from here on referred 
to as chiral 2DEGs (C2DEGs)^. C2DEG models pro- 
vide an accurate description of the low-energy properties 
of few-layer graphene systems with a variety of different 
stacking arrangements consistent with symmetries. 

Properties of quasiparticle excitations in C2DEGs are 
determined by their chirality index J. The quasiparticle 
dispersion is given by ej ~ |p|"' (where p is momentum 
quantum number measured about some special points 
in the Brillouin zone) and quasiparticles exhibit a Berry 
phase of Jtt. These properties lead to an unusual Landau 
quantization. Due to this, integer quantum Hall effects 
(IQHEs) in C2DEGs are remarkably different from that 
of semiconducting 2DEGs. Examples of this are quantum 
Hall effects in single layer, bilayer and trilayer graphene 
which exhibit unusual IQHEs described by Dirac con- 
tinuum models. At low-energies these systems represent 
the J = 1, J = 2 and J = 3 instances of the C2DEGsS 
model respectively. Another unique property of C2DEGs 
is the presence of a zero-energy LL with degenerate LL 



orbitals n = O,-- - , J ~ 1. Speaking loosely, quantum 
states corresponding to cyclotron orbits with different ra- 
dius, which would have different energies in an ordinary 
two-dimensional electron gas, are degenerate. This de- 
generacy is, of course, on top of the normal Landau level 
one- state-per- flux- quantum degeneracy and the four-fold 
degeneracy already present due to spin and valley degrees 
of freedom. The zero-energy LL is thus 4 J-fold degener- 
ate. 

Naively it can be anticipated that Coulomb interac- 
tions will lift the 4 J-fold degeneracy of the A^ = LL by 
producing spontaneously broken-symmetry ground states 
with spin, valley pseudospin and LL orbital pseudospin 
polarizations. As a consequence, one would expect quan- 
tum Hall plateaus at all intermediate integer values of 
the filling factors from i/ = — 2J+1 to i' = 2 J—"—. It 
turns out that this expectation only holds for J < 2 
C2DEGs, whereas for J > 2 insulating states in the top- 
most LL appear generically, which leads to a disappear- 
ance of certain Hall plateaual? from the above sequence 
and a quantization of the Hall conductivity cr^y at a value 
corresponding to the adjacent interaction driven integer 
quantum Hall plateau (for example, axy = —6e^/h when 
ly — —5). Thus, the presence of wave functions with 
different spatial structures appearing at the same energy 
has important consequences on the interaction driven QH 
states in ABC-trilayer graphene (i.e. C2DEG systems 
with chirality index J > 2). In particular due to the LL 
orbital degeneracy, spin and valley polarized non-uniform 
states are allowed due to charge modulation in the orbital 
subspace. As we show these states only appear for J > 2 
C2DEGs, and are the main focus of this paper. 

The existence of these additional plateaus from sponta- 
neously broken-symmetry ground state has already been 
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confirmed in suspended bilayer graphene samples and 
bilayer graphene on Si02/Si substrates^'^ as well as in 
graphene trilayersi^. By applying an electric potential 
difi^erence (or bias) between the outermost layers 
(where the low-energy sites reside) of a graphene bi- 
layer, one can control the population of electrons in each 
layer and open a gap All [As, B) between two adja- 
cent orbital state which is a function of both A b and the 
quantizing magnetic field B. The phase diagram of the 
C2DEG in the v — All space has been studied for bi- 
layer graphene and is very rich. States with spin and/or 
layer and/or orbital polarizations are possibleiSii^ii^^Ji. 
At V = —1,4, the bias drives a series a transitions from 
an homogeneous state with finite orbital pseudospin to 
charge-density-wave states and crystal states where the 
orbital pseudospin rotates in spacei^»i2,. Interestingly, 
a similar sequence of transitions has been observed in a 
thin film of the helical magnet Feo.5Coo.5Si using Lorentz 
transmission electron microscopyi^^— . In that system, 
the phase transitions are induced by a transverse mag- 
netic field. In the graphene bilayer, they are induced by 
Ab potential. The Hamiltonian of both systems are sim- 
ilar. In particular, they contain a Dzyaloshinsky-Moriya 
(DM) interactional which, in the bilayer, is entirely due 
to exchange interactions instead of from spin-orbit cou- 
pling. 

In this paper, we study the quantum Hall ferromag- 
netic states of the C2DEG in trilayer graphene's zeroth 
LL. Assuming the layer and spin degrees of freedom are 
inactive^^, we only consider the orbital pseudospin at 
integer filling. This situation occurs at filling factors 
1/ = —5,-4 and = 4, 5 for finite value of lA^^I or at 
other filling factors in some specific range of All- The 
phase diagram consists in a uniform phase which appears 
at large absolute values of All where electrons occupy 
the lowest energy levels and there is no orbital coher- 
ence. Upon decreasing the value of |All|, the uniform 
state shows an instability in the pseudospin wave mode at 
finite wave vectors q . This instability indicates, in prin- 
ciple, a transition to a unidirectional charge-density-wave 
state. We find, however, that this transition is preempted 
by a first-order transition to a crystal state with a vor- 
tex pseudospin texture. One remarkable result is that 
this crystal state in the trilayer exists at both positive 
and negative values of All in contrast with the bilayer 
case where it is found for A^l < only. The pseudospin 
texture of the crystal is more complex than in bilayer 
graphene since an electronic state in this system is de- 
scribed not by a CP^ but by a CP^ spinor. This spinor 
can be decomposed into three distinct pseudospin vec- 
tors. We show that our crystal has a vortex texture in 
each of these pseudospins but with different vorticities. 
We discuss some properties of the crystal and uniform 
phases in this article. In particular, we show that these 
states can be distinguished by their electromagnetic ab- 
sorption spectrum. 

Our paper is organized in the following way. In Sec. 
II, we review the effective two-band model Hamiltonian 



for ABC-trilayer graphene. In Sec. HI, we describe the 
Hartree-Fock and generalized random-phase formalism 
(GRPA) used to compute the energy and pseudospin tex- 
tures of the different phases as well as their collective 
excitations. In Sec. IV, we explain our pseudospin repre- 
sentation for the various phases, define the electric dipole 
density and the formalism used to compute the electro- 
magnetic absorption. The phase diagram of the C2DEG 
is presented in Sec. V. In Sec. VI we show that the uni- 
form phase is unstable at a finite wave vector in some 
range of bias where the crystal state described in Sec. 
VII emerges. We discuss the properties of this crystal 
state in Sec. VIII and conclude in Sec. IX. Appendix A 
lists the values of the Coulomb exchange interactions at 
zero wave vector. Appendix B summarizes the Hartree- 
Fock and GRPA equations for the single and two-particle 
Green's functions and Appendix C gives the generators 
of the group SU(3). 



II. ABC-TRILAYER GRAPHENE AS A C2DEG 

In this section, we describe the low-energy effective 
Hamiltonian for the C2DEG in ABC-stacked trilayer 
graphene. Before we look at the specific details of ABC- 
stacked trilayer graphene lets us review the properties 
of C2DEGs^. The low-energy effective Hamiltonian of 
a chiral two-dimensional electron gas (C2DEG) can be 
written^ as 



[cos {J6) + sin {J9) ay] , (1) 



where axjCTy are Pauli matrices, p is the momentum of 
the electron, 9 its angle with the x axis and J is the chi- 
rality index. The parameter pc = 71 /vq where 71 is the 
interlayer hopping energy between the two high-energy 
sites in adjacent planes and vp = 3co7o/2fi is the Fermi 
velocity with cq the separation between carbon atoms in 
a plane and 70 the intralayer hopping energy between 
two neighboring carbon atoms. This Hamiltonian oper- 
ates in the space of a two-component wave functions 
describing electronic amplitudes on the two low-energy 
sites A and B. In the valley K = (-2/3, 0) 27r/a, C = -1 
and = (A) , i/j (B)) whereas in the valley K' = 
(2/3, 0) 27r/a, ( = +1 and *+ = (V' (B) , ^ (A)) . (For an 
introduction to the electronic properties of C2DEG's, see 
Ref. 8). 

When a transverse magnetic field is applied to a 
C2DEG, the kinetic energy of the electrons is quantized 
into Landau levels with energies 



En,^ = sgn (N) 71^ 



hi 



j-i 



11(1^1-*)' (2) 



where sgn is the signum function and i — yJhc/eB is 
the magnetic length, N = 0,±1,±2, ... is the LL index 
The LLs are 4-fold degenerate when counting valley and 



3 



spin degrees of freedom. The = LL is special since, 
including the valley and spin degrees of freedom, it is 4 J- 
fold degenerate. The extra degeneracy comes from the 
fact that the eigenspinors in = which have the form 
{( hn,x (r) ) , n = 0, 1, J — l} are degenerate (here, 
hn,x (r) represents the Landau-gauge wave functions of 
conventional 2DEGs given in Eq. (|6|) below). We refer to 
the index n as the orbital quantum number. The pres- 
ence of wave functions with spatial structures that ap- 
pear at different energies in the ordinary non-relativistic 
2DEG model in the same degenerate manifold can create 
some terminological confusion. We will refer to the wave 
functions hn as Landau level n orbitals and (as already 
anticipated) use upper case letter N to distinguish lev- 
els with different Landau-quantized band energies in the 
C2DEG model in a magnetic field. 

The lattice structure of an ABC-stacked trilayer 
graphene (rhombohedral stacking) is shown in Fig. [T] 
Each layer has a honeycomb lattice of carbon atoms. The 
underlying Bravais lattice is a triangular lattice with a 
basis of two atoms denoted by Am and Bm where m is the 
layer index. The triangular lattice constant is ag = -s/Scq 
where cq = 1.42 A is the distance between two neighbor- 
ing carbon atoms. The Brillouin zone of the reciprocal 
lattice has two non-equivalent K points that we take as 
K± — ± (2/3, 0) 27r/ao as indicated in the inset of Fig. 
[TJ Each adjacent layer pair forms an AB-stacked bilayer 
with the upper B sublattice directly on top of the lower A 
sublattice. The upper A sublattice is above the center of 
a hexagonal plaquette of the layer below. Two adjacent 
layers are separated by a distance d = 3.35 A. 

The band structure of the ABC-stacked trilayer 
graphene has been studied in Refs l22ll23l . Near the val- 
leys K±, it consists in three valence and three conduc- 
tion bands as shown in Fig. (3] In the simplest model 
where only the nearest-neighbor intralayer 70 ~ 3.16 eV 
and interlayer hopping 71 w 0.502 eV are considered, the 
degenerate bands in the middle of Fig. [2] have a cubic 
dispersion. For undoped ABC-trilayer graphene, the va- 
lence bands are completely filled and the Fermi level lies 
at ii^ = 0. The high-energy bands are separated by a gap 
71 from the low-energy bands as shown in Fig. [2] The 
low-energy bands touch at the K± points while the other 
four bands cross at the energies E — ±71 above (below) . 

To study the low-energy behavior of the electrons, we 
use an effective two-band model which results from per- 
turbation theory in VQ/jf. This effective model can be 
derived for ABC-stacked trilayer^^'^- starting from a cou- 
pled Dirac model, consistent with the stacking arrange- 
ment. In the basis of the low-energy sites {Ai , B^) for the 
valley K_|_ and (B^jAi) for the valley K_, the resulting 
Hamiltonian, in the case where a perpendicular magnetic 
field is applied, is given by 

- ^ 3 ^ L (3) 




FIG. 1: (Color online)Lattice structure of an ABC-stacked 
graphene trilayer. The two non-equivalent sites of the honey- 
comb lattice in each plane are indicated by Am and Bm , where 
m is the layer index. The two basis vectors of the underlying 
hexagonal Bravais lattice are ai,a2.The Brillouin zone of the 
hexagonal lattice with the two non-equivalent points K± is 
drawn in the top right corner of the figure. 
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FIG. 2: Band structure of ABC-stacked trilayer graphene 
obtained from the tight-binding Hamiltonian by keeping the 
hopping parameters 70 and 71 only. 
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where 



(4) 



and /3 = vo/ji.ln Eq. (jS]), a, are the ladder operators 
for the Landau levels. The valley index is ^ = ± and 
we have defined Vi — y/3/2ao^i/i. In deriving, we 
have taken into account a perpendicular electric field that 
forces a potential difference (or bias) between the 
outermost layers. In a magnetic field, the low-energy 
bands of the full model are replaced by a set of Landau 
levels with energies, in the absence of bias, given by Eq. 
([2]) with J = 2. These energies are independent of the 
guiding-center coordinate X so that each level has the 
usual degeneracy = S'/27^£^ where S is the area of 
the 2DEG. The three eigenspinor for the orbital states in 
= and in the Landau gauge A = (0, Bx, 0) are given 

by 





^2,x (r) 





hi,x (r) 





ho,x (r) 



wher( 



.24 



(r) 



(5) 



(6) 



are the wave functions in the Landau gauge with ifn (x) 
the wave functions of the one-dimensional harmonic os- 
cillator. 

A finite bias lifts this orbital degeneracy. The energies 
are proportional to and are given as 



where 



,A, 



(7) 



(8) 



The correction All is small compared to the bias. In- 
deed, if we use the values of the tight-binding parameters 
given in Ref. we find — 5.49 x lO^^B where B 
is the magnetic field in Tesla. One remarkable aspect of 
the energies of the orbital states is that the ordering of 
the energy levels in A = is different in the two valleys. 
This property of the energy spectrum has profound con- 
sequences on the phase diagram of the C2DEG in trilayer 
graphene, as we will show in this paper. 

Other remote inter-layer hopping parameters that have 
hitherto been neglected in our analysis can modify the 
energy spectrum. The hopping term 74 (which couples 
the low- and high-energy sites located on different lay- 



ers) adds a finite correction —2n(3i/4 to E9 



which 



is independent of the valley index and bias and scales 
linearly with the magnetic field. This correction lifts the 
degeneracy of the orbital states even at zero bias. We 
can include it in All by redefining 



,A, 



2^/31^4. 



(9) 



In this paper, we take A^i (not A^) as the parameter 
that we vary to study the phase diagram of the C2DEG. 
Clearly, All can be tuned by changing the bias or the 
magnetic field. It can have both positive and negative 
values. The maximal value of |All| must be such that we 
stay within the limit of validity of the two-band model. 
This can be checked by comparing the band structure 
of the two-band model (Eqs. dTHHl)) with that given 
by the full (six band) model with all hopping terms in- 
cluded. We have done this comparison and will report it 
elsewhere^^. Our conclusion is that there exists a range 
of bias where the two-band model is well-justified. This 
range increases with increasing magnetic field. 



III. ORDER PARAMETERS AND COLLECTIVE 
EXCITATIONS 



In the present work, we study the phases of the C2DEG 
when the Fermi level, filling factor and bias are such that 
the trilayer can reasonably be described by a three-level 
system in the N = LL with level energies given by 



LL- 



(10) 



To do this, we must consider the valley and spin degrees 
of freedom to be frozen. This can occur, for example, 
at filling factors = —5, —4 or at = 4, 5 when the 
lower levels are fully filled and can be considered as in- 
ert. Also, when Coulomb interaction is included, the Zee- 
man gap is exchange-enhanced and, for the filling factors 
just mentioned, the ground states were shown to be spin 
polarized^^. Finally, interlayer coherence occurs only for 
very small bias |All| ^ 0.001 (? j (we checked this 
numerically) so that, unless All is close to zero, layer 
polarization can safely be assumed. 

We denote by Vn the filling factor of the orbital level n 
and by v the filling factor of the three-level system. Our 
aim is to study the phase diagram of the C2DEG when 
V — (since 1/ = 3 is trivial) as All is varied. 

To study the phase diagram, including both homoge- 
neous and modulated states, we define the operators 



(11) 



X.X' 



t JC 
'^^n,X^'n-' ,X'Ox,X'+qye^i 



where A^ = Sj^-Kp' is the Landau level degeneracy. 
The Hartree-Fock Hamiltonian can then be written as 



Hhf = A^^J!p„,„ (0) (12) 

(q) (p ni,n2 V q))pn3,n4 (q) 

q 

'-^V ^ -'^ni,n4,n3,n2 (q) {Pni;n2 (~q)) Pn3,n4 (q) J 
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where repeated indices are summed over. In deriving Eq. 
(|12p . we have taken into account a neutrahzing positive 
background so that the q = contribution is absent from 
the Hartree term. This is indicated by a bar over the 
summation. 

The Hartree and Fock interactions are defined by 



X 



ni,"2.n3,"4 



(q) = 



H. 



ni ,712 i''^'' 



The Coulomb energy e^/nl = 56.2\/B meV with 
B in Tesla and k = 1. The Fock interactions 



X 



ni,n2,"3,n4 



(q = 0) are Hsted in Appendix A. 



by 



The form factors which appear in H and X are given 



then the order parameters of the coherent phases are sim- 
ply 



(P«i,«2 (q)) = Gn^^ni (q,r = ) 



(20) 



The equation of motion for the Green's function in the 
Hartree- Fock approximation is given in Appendix B. This 
equation leads to the sum rule (at T = K) 

EEi(^— = (o)>. (21) 

q "2 

By definition, we also have 

{Pn.n (0)) = Vr,. (22) 

To study the collective excitations, we compute the 
two-particle Green's function 



Kn,,n2 (q) = 



with 



Fn,n' (q) 



Fni,n2 (q) if ni > n2 
[Fn2.n1 (-q)]* if ni < n2 



In']. ( {qy + iqx)i 
,1 



V2 



(15) 



(16) 



rn-n' 



They capture the character of the different orbital states. 
Finally, the Hartree-Fock energy per electron is given 

by 



^ = IeI (p„,„ (0)) 



(17) 



2V 



E^"i 



"2, "3, "4 V^J \rni,n2 \ 

n4,ri3,rt2 V q) (P ni ;n2 (-q)) (p«3,n4 (q)) , 



where is the number of electrons in the C2DEG. 

The order parameters of the orbital phases are ob- 
tained from the single-particle Matsubara Green's func- 
tion 

G„, (X, X', r) = - (t.c„„x (t) cl^^, (0)) , (18) 

where Tr is the imaginary time ordering operator and 
c]j ^ creates an electron in orbital n with guiding-center 
X.' 

If we define the Fourier transform of the single-particle 
Green's function as 



G 



..«2(q,r) = ^ E 



(19) 



X,X' 



(X,X',t) 



Xni,n2.n3.n4 

(q,q';r) (23) 

= -N^ (7VPni,«2 (q,'r) Pn3,n4 (^q'; 0)) 

ni ,7i2 (q)) {pns.m (-q')> 

in the generalized random-phase approximation— 
(GRPA). The resulting set of equations is 
given in Appendix B. The collective exci- 
tations are given by the poles of the re- 
tarded Green's function Xrn'!n2,n3,"4 (qi qi <^) = 
Xni.n2,n3.n4, (q, q; i^n uj + iS) .To derive the dis- 
persion relations, we follow these poles as the wave 
vector q is varied in the Brillouin zone. 



IV. DESCRIPTION OF THE 
ORBITAL-COHERENT PHASES 

A. Pseudospin representation 

An electronic state in our three-state model can be 
created by the spinor field 



$t (r) = I vl/t (r) 
vft (r) 



(24) 



where ^-t (r) = /i^ (r) cl ,,. 

Using the eight infinitesimal generators Ta (with a = 
1,2,..., 8) of SU(3) (see Appendix C), we can define the 
eight real fields 



Fa (r) = $t (r) (r) 
with the Fourier transforms 



Fa (q) = J dre-^'i--$t (r) T,<i> (r) , 



(25) 



(26) 
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With a = 1, we get 

^1 (q) = \ [N^Ko,i (-q) po.i (q) + A^0^i,o (-q) Pi,o (q)] • 

(27) 

We describe the three-level system by three 
pseudospins^^ 1/2. We associate spin up (down) 
with level n — to get the spin system. We take 
(ij) = (0, 1) ; (1,2) ; (0,2) . We suppress the orbital- 
dependent part of the form factor and keep only the 

factor /3q = e ~ in Kij (q) for aU i,j. Eq. (|27|) . for 
example, becomes 



i^i(q)= V^i) (q). 



(28) 



where pi* '"'^ = ^ iPi.j + Pj.i) ■ T^^^ other components are 
listed in Appendix C. We also add to these fields the 
"densities" 



Pt (q) = PqPia (q) , 



(29) 



with i = 0,1,2. 

The eight real fields Fi (r) (the Fourier transforms of 
Fi (q)) provide a complete description of each phase stud- 
ied in this paper. 



B. Dipole density 

The total electronic density is given by 

n{r)^Y.^Ur)^,{r). (30) 

Its Fourier transform is 

2 

n (q) ^^'1 ("^) P^'i ■ (31) 

An external electric field, Eext (r) — — V0ea;t (r) , couples 
to the density through a term 



in the Hamiltonian where 



(r) = ^^JZ"^-* 



(32) 



(33) 



Using the definition of the form factors Ktj (q) given 
in Eq. P3|) . the coupling iJe^t can be written as 

H,,t = J drpTOT (r) (b{r)- J dr (d (r) • E (r)) , (34) 



where 



PTOT (q) = -e- J2 N^e-^^'^^'ep.,, (q) , (35) 



and we can define the dipole operators 

4(q) = -7(q)U°'^)(q) + «(q)p^2)(q) , (36) 



dyiq) = 7(q) 



pf^) (q)-Fa(q)p(i'2) (q) 



^/2 (1 - g2^V4) and 7 (q) 



(37) 



with a (q) = 
\/2iV^efe-9'^V4. 

In the phases studied in this paper, {proT (r)) is al- 
ways uniform i.e. X)i=o (Pi^i (q)) ~ ^'^q.o for = 1,2. It 
follows that we can ignore the first term in H^xt in Eq. 
([55]) and the coupling with the external electric field is 
simply 



ffe,t = -y"dr(d(r)-E(r)) 



(38) 



where d (r) can be interpreted as a density of electric 
dipoles^*. 

In the absence of Coulomb interaction, the time vari- 
ation of the total dipole moment is given by 



..ld(0) 



-[<^,d(0)] 
= iALLZ X d (0) , 

where 

H^HF - [AllPi,i (0) + 2AllP2.2 (0)] 
so that the dipoles oscillate at the frequency 

ujdip = All/H. 
We can define a dipolar current density by 



J*P = -d (0) . 



C. Optical absorption 

The total current in the 2DEG is given by 
1 



(39) 



(40) 
(41) 

(42) 



J = dr 



(*t (r)jvI/(r)) + (jvI/(r))tM/(r) 



where 



(43) 



(44) 



with A*^ the vector potential of the external electromag- 
netic field and H'^ is the Hamiltonian of the two-band 
model. We find 



Jx = -4S 



(0) + V2p(i-2) (0)' 



i=0 



Jy = -4S Uo.i) (0) -f V2p(i'^) (0) 



(45) 
(46) 
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with the constant 



1 



A 



(47) 



The total current^^ given by Eqs. (|45ll46p is nothing but 
the dipolar current defined in Eq. (|^^ above i.e. J = 

J dip 

The optical absorption per unit surface from an elec- 
tromagnetic wave E ~ EoSae^'^* (with a — x,y and we 
define a = x if a — y and a = y if a — x) is obtained 
from the retarded current-current response function 



n 



Lo -\- i5 



(48) 





2 

j Im 


'Xoi.-a (0,t^)" 






uj + i5 



where the retarded current response function is obtained 
from the time-ordered two-particle Green's function 



with 



and 



XJ.J, (r) = - (TJ, (r) Jp (0)) , (49) 



Xcp (0, r) = - (Tp„ (0, r) (0, 0)) (50) 



p„(0,r) = p(O-i)(0,r)-f ^/2p(i-2)(0,r). (51) 



V. PHASE DIAGRAM OF THE C2DEG 

The ground state in the ABC-trilayer graphene can be 
classified in terms of translationally invariant (uniform) 
states or non-translationally invariant (non-uniform) 
states. One important distinction to note is that any 
non-uniform state in graphene would otherwise be either 
spin or valley density waves. However, for spin and valley 
polarized J > 2 C2DEGs, non-uniform states are gener- 
ically allowed due to charge modulation in the orbital 
subspace. Both states have different experimental signa- 
tures: uniform states exhibit Hall conductivity whereas 
if the topmost occupied LL has a crystal-like state or 
unidirectional charge-density-wave, it will either be in- 
sulating or exhibit anisotropic conductivity. In the case 
of crystal-like states discussed here, which will likely be 
pinned by disorder, the Hall conductivity a^y will be at 
a value corresponding to the adjacent interaction driven 
integer quantum Hall plateau as we will discuss. 

In our study of the phase diagram, we consider the 
following states: 

1. A coherent uniform state (CUP). In this state, the 
only allowed order parameters are (yOn,m (q = 0)) . 
The state of each electron is described by the CP^ 
spinor (00,01,02) (where a^'s are complex num- 
bers satisiying l'^"!^ = 1) so that an electron 



at guiding-center X is in a linear combination of 
the three orbital states. This combination is the 
same for all electrons. The CUP ground state is 
written as 



i*)ct/p=n 



X 



X 



|0) 



which gives 



{Pn,m (q = 0)) = a* a„ 



(52) 



(53) 



Due to particle-hole symmetry the ground state for 

V — 2 can be described as a filled level of holes on 
a vacuum state consisting of the three levels filled 
with electrons. The CP^ spinor (00,01,02) also ap- 
plies to an hole state if — >■ h\-^ (where h^^^ 
is a hole creation operator). The CUP phase is 
possible at negative bias because the system can 
then reduce its kinetic energy by populating the 
levels n = 1,2 that are below level n = in en- 
ergy. This, however, increases the exchange en- 
ergy because the Coulomb exchange terms satisfy 
^0,0,0,0 (0) > ^1,1,1,1 (0) > ^2,2,2,2 (0) (see Ap- 
pendix A). In consequence, there is an optimal pop- 
ulation of the levels that minimizes the total energy. 

2. An incoherent uniform phase (lUP). In this 
case, the only allowed order parameters are 
{Pn,n (q = 0)) and the first { u = 1) oi first two 
{v — 2) lowest-lying orbital states are fully filled 
so that all coherences a'^a^ {n ^ m) are zero. For 

V — 1, we have in this limit oo = 1 when A^^ > 
and 02 = 1 when A^^ < 0. (It is just the opposite 
for the hole spinor when v = 2.) 

3. A coherent charge-density- wave phase (CCDWP). 
This state is modulated in one direction only 
and the allowed order parameters are {pn,m (pqo)) 
where p — 0, ±1, ±2, ... and qo is the wave vector of 
the CCDWS. The ground state is simply obtained 
by letting o„ o„ (A") in Eq. ((5^ . 

4. A coherent crystal phase (CCP). In this non- 
uniform state, all order parameters {(p„_m(G))} 
are allowed where {G} are the reciprocal lattice 
vectors of the Bravais lattice of the crystal. We 
have considered a triangular and a square lattice 
with one electron per unit cell. 

Note that in all of these states, the "density" {p (r)) = 
X^n {Pn,n (r)) IS a constant in space. The real density 
{n (r)) as defined in Eq. (pO)) is modulated in space in 
the CCDWP and CCP. 

Our Hartree-Fock numerical calculations for the phase 
diagram of the C2DEG is shown in Table 1. The CUP 
and CCDWP are never the ground state. The ground 
state is an lUP for large value of lA^^I and a CCP (with 
a triangular lattice) in between. The transitions from the 
lUP to the crystal state are first order. 
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Phase 


All 


lUP in n = 2 
CCP 

lUP in n = 


All < -0.31(-0.34) e^/^l 

-0.31(-0.34) e^/K£ < All < 0.09(0.03) e^/ni 

All > 0.09(0.03) e'^/nl 



TABLE I: Phase diagram including the uniform and crystal 
phases for v — 1. The numbers in parenthesis are for v — 2. 



We discuss the crystal state in more details in the next 
section. We point out here, however, that the presence 
in the phase diagram of a crystal state at positive value 
of is something specific to trilayer graphene or other 
C2DEGS with J > 2 as discussed in Sec. IVira below. 



VI. PROPERTIES OF THE UNIFORM PHASES 

Even though the CUP is not the ground state in tri- 
layer graphene, we would like to give a brief description 
some of its properties. In particular, it shows an instabil- 
ity to a CCDW state that is, in the case of the trilayer, 
preempted by the crystal phase. Comparing the ener- 
gies of the lUP and CUP in the HFA, we find that the 
CUP has lower energy than the lUP in the range of bias 



All G 



(*) 

LL 



where A^*l = -0.29e'^/Ki for i> = 1 and 



-0.24e^/Ke for v 



A 

The ground-state energy per electron is an uniform 
phase is given by 



— — = - [All (Pi.i) 



-'^o, 0,0,0 (po,o) 



1 

- [-'^o, 1,1,0 (Po,o) (Pl,l 

V 



f2Aii(p2,2)] (54) 

-'^l, 1,1,1 (Pl,l)^ + -^^2,2,2,2 {P2,2Y 
) + XQ,2,2fi (P0,0> (P2,2)] 



~ [-''^1,2,2,1 (Pl,l) 



(P2,2)] 



0,0,1,1 



l(PO,l 



V 

— [-''^0,1,2,1 (Po,l) 
V 



{P2 



-^^0,0,2,2 |(P0,2)| 
> + -^^1,0,1,2 (Pl,2) 



f -^^1,1,2,2 |(Pl,2) 

(Pi,o)] 



where all (pij) ' s and the interactions Xi jj^ i are evalu- 
ated at q = 0. Note that if we make a local gauge trans- 
formation of the CP^ spinor i.e. a; {X) — )• aie^^'-^\ the 
energy Ecup is unchanged. This U{1) gauge invariance 
is necessary to define the CP^ spinor-30. Fig. |3]shows the 
occupation of the levels in both phases (lUP and CUP) 
for filling factors = 1,2. 

One possible parametrization of a CP^ spinor is given 
bj^ 



cosf 
e*" sin 9 cos ip 



(55) 



AJ{e'/K\) 




AJ{eVK\) 




AJie^M) 

FIG. 3: (Color online) Occupation of the orbital levels in 
the incoherent {All > 0) and coherent {All < 0) uniform 
phases for _B = 10 T as a function of the bare gap A_LL.(a) 
1/ = 1 and (b) 5 = 2. 



sm f sm If 



The first four lines on the right hand side of Eq. ([54]) are 
independent of the angles a and /3 while the two terms 
in the last line depend on cos {a — /3) . We have verified 
numerically that the ground state energy is minimized 
when a — P and that it is independent of the choice of 
a. The CUP thus has a broken U(l) symmetry and sup- 
ports a Goldstone mode. This is confirmed by our GRPA 
calculation which also shows that the dispersion of this 
mode is highly anisotropic. The CUP in a graphene bi- 
layer has similar properties as reported in Ref. [VQ where 
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the origin of the anisotropy is discussed. 

By contrast, the lowest-energy mode in the lUP is 
gapped and has an isotropic dispersion. This gap can be 
calculated analytically from the GRPA equations. We 
find 



0) 



(56) 



15 
64' 



Ye! 

\^ Kt 



if J7= 1, 

if 1, 



All > 
All <0 



uiiup (q ^ 
' All, 
-All 

All, if 1^=2, All > 

-ALL-^,^{^),if^^2,ALL<0 

The critical bias Aj^£ for the transition from the lUP 
to the CUP is given by the condition ujiup (q = 0) = 0. 
The frequency loiup (q = 0) is however positive in the 
region where the lUP is the ground state according to 
Table 1 . This frequency is measurable in electromagnetic 
absorption experiments. We come back to this point in 
Sec. VIII. 

The dispersion of the lowest-energy mode in the lUP 
and CUP becomes unstable at a finite value of q in some 



range of bias All ^ 



^(1) ^(2) 



.Fig. m shows the sit- 



LL' ^LL 

nation for = 1. In this case, the lUP is unstable for 
(j^ sa 2 at A^^l = 0.016 e^/k^ while the CUP is unstable 

for qyi « -2 and qyi « -3 at A^^| = -0.25 e'^/ni. The 
direction in q— space of the instability is related to the 
orientation of the electric dipoles present in the CUP. The 
instability in to (q) occurs in the direction q = z x d (r) 
as in the bilayer case^^. The dispersions at these two 
biases are plotted in Fig. [SJ For filling factor ? = 2, it 
is the lUP at negative All that becomes unstable for 
All > A^"'^ — —0.27 e^/ni and the instability persists 



well into the CUP until A 



(2) 
LL 
(2) 



'0.036 c'^/kI. The sys- 



— ■ — \ 


L=0.016 


— • — \ 


L=-0.25 



.2. 




tem is stable for All > A}^^. 



FIG. 5: (Color online) Dispersion relation of the collec- 
tive modes in the uniform phases near the instability points 
A'-^I = -0.25 e'^/Ke (CUP) and A^^^ = 0.016 e'^/Ki (lUP) 
defined in Fig. 4. 



of some kind of charge-density-wave state. Some of us 
discussed this instability in ReffT^ For negative value of 
All, this is not surprising. The same situation occurs in 
bilayer graphene^^^. The instability for positive value 
of All is unexpected, however, and does not occur in 
bilayer graphene^^. 



VII. CRYSTAL STATE WITH ORBITAL 
COHERENCE 



lUP 




CUP 


lU 


P 






unstab 


le 













A A 



A 



(2) 
LL 



FIG. 4: (Color online) Phase diagram ai v = 1. The in- 
coherent uniform phase (lUP) occurs for All > and for 
All < ^LL- between these two biases, the system is in 
a coherent uniform phase (CUP) in the HFA. In the GRPA, 
the collective excitations show that the uniform phases are 
unstable between A^^^ A^^^- 

The instability at a finite wave vector suggests that 
both uniform phases are unstable towards the formation 



In the HFA, we find that the instability towards the 
CDW state is preempted by the formation of a coherent 
crystal phase (CCP) with a triangular lattice. This elec- 
tron crystal has one electron (i/ = 1) or one hole = 2) 

per site so that the lattice constant is cq /£ = In j ^/^J^ 
in both cases. The real density n (r) is modulated in 
space but not the "density" (p (r)) = ^„ (pn,n (r)) which 
is a constant. At integer filling, this is possible because 
in our system, the electrons can be distributed in more 
than one orbital states. Every electron in the CCP is in 
a linear superposition of the three orbital states as in the 
CUP but also in a superposition of guiding center states 
X. All order parameters (pn,,n (G)) are finite where {G} 
are the reciprocal lattice vectors of the crystal. We choose 
I G I big enough in the numerical calculation to insure that 
the crystal energy converges to the required accuracy. We 
show in Fig. [S] the filling factors z/„ of the orbital states 
as well as the cohesive energy of the crystal. This energy 
is quite large, of the order of 0.04 / kI w 7 meVw 80 K 
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at B = 10 T and for k = 1. Hence, this state should be 
quite robust against thermal fluctuations and disorder. 



0.010 




ALL/(e^/Kl) 



FIG. 6: (Color online) Filling factors (left y axis) of the orbital 
levels n = (squares), n — 1 (delta), and n — 2 (nabla) in the 
cystal states for v = 1 and B = 10 T. The difference in energy 
between the crystal and the uniform states is shown by the 
green curve with the right triangles (right y axis). The dashed 
vertical line indicates the point where the uniform state is 
lower in energy than the crystal state for All positive. 

We show in Fig. [7] the electronic density n (r) and 
the full SU(3) pseudospin representation of the crystal 
state for i/ = 1 at S = 10 T. We have chosen All = 
in this figure, but the textures do not change much as 
All is varied. We remark that the fields Fa (r) defined 
in Eq. (j28p represent pseudospin densities. These are 
not bounded. In particular, they are not normalized in 
Fig. [71 We find that for the three pseudospin fields, 

(4*'^') (r)) » ^f^'''> (r)) , (r)) so that the pseu- 

dospin are almost completely polarized in the direction 
of the z axis or opposite to it. The in-plane compo- 
nent of the pseudospin vectors (and so the interorbital 
coherence) is small. There is hovirever a small clockvifise 
rotation of the in-plane component of the pseudospins in 
the x — y plane around each lattice site. For (i, j) — (0, 1) 
and (1, 2) the rotation of the pseudospins is of 27r while 
for (0, 2) the pseudospins rotate by An. 

We conjecture that the interesting pseudospin texture 
around each lattice site can be assigned a topological 
charge. The minimal CP^ sigma model is known to sup- 
port skyrmion solutions^. In our case, where multiple 
orbitals are considered, the energy of long-wavelength 
deformations contains much more terms than the min- 
imal CP^ sigma model. In particular, it involves multi- 
ple pseudospin stuffiness. Nevertheless, we believe that 



finite energy excitations should fall into different topo- 
logical sectors even in this case. More work is needed, 
however, to confirm our conjecture. 

A more physical quantity to represent graphically is 
the dipole density defined in Eqs. (136137^ . It is shown 
in Fig. [8] for filling factors = 1,2. The vector field 
of the dipole density has also a vortex structure around 
each lattice site. The dipole vectors rotate by 27r in both 
cases. 



A. Absence of the crystal state in bilayer graphene 

As shown in Fig. [3J the effective interactions 
Hni,n2,n3,ni (q) " -'^«i,n4,«3,n2 (q) which appears in the 
Hartree-Fock energy favor the formation of a CCDWP or 
CCP state because most of them take their minimal value 
at a finite wave vector. Because the filling factor is an 
integer, the system must put some of the electrons in the 
higher-energy orbital states in order to produce a density 
modulation. These modulations increases the Hartree en- 
ergy of the system and, in the case where All > 0, the 
occupation of the higher-energy levels increases the bias 
energy. Nevertheless, the crystal state is favored in the 
trilayer because these costs are more than compensated 
by the gain in exchange energy as shown in Fig. 1101 

In the Bernal-stacked bilayer, the Hartree-Fock equa- 
tion admits a crystal solution for All > but its total 
energy is greater than that of the lUP (see [TU] (b)). We 
believe that the crystal energy in the bilayer case is not 
optimal because of the following reason. The effective in- 
teractions are not monotonous as shown in Fig. [9l At the 
value qo where the energy is minimal, the particular ef- 
fective interaction i?„,„,„,„ (qo) - Xn^n,n,n (qo) is smaller 
for n — 0. At larger value of q, it is the opposite i.e. 
the interaction is smaller for n = 2. (The ordering in en- 
ergy of the interaction i?ni,n2,n3,n4 (q) - ^ni 
depends very much on the indices ni, ...,^4). In a crys- 
tal, {pni,n2 (G)) is non zero for an infinite set of val- 
ues of G. These order parameters are constrained by 
the sum rules of Eq. (PT|) and also by the condition 
J2n {Pn,n (r)) = J' . In thc bilayer case, this condition 
imposes (po,o (G)) — — {pis (G)) for G ^0. This severe 
constraint does not allow the crystal state to take full ad- 
vantage of the non-monotonous behavior of the effective 
interaction in distributing its weight amongst the order 
parameters with different G's. By contrast, the trilayer's 
constraint (po,o (G)) -I- (G)) + (p2,2 (G)) = is much 
less severe. This point has been discussed on symmetry 
grounds in Ref. |l3- The effective theory describing the 
CCP to lUP phase transition requires the presence of 
third order term which, based on symmetry arguments, 
vanish for any two-level system such as bilayer graphene. 
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FIG. 8: (Color online) Electronic density and dipole field in 
the crystal state for (a) = 1 and (b) v = 2. Parameters are 
All = and B = 10 T. 



VIII. PROPERTIES OF THE CRYSTAL STATE 

In this section, we look at some of the properties of the 
coherent crystal state in more details. 



A. Density of states 



FIG. 7: (Color online) Pseudospin texture in the (a) (0, 1); 
(b) (0,2); (c) (1,2) pseudospin systems for the crystal state 
at 1/ = 1, All = and i? = 10 T. The total electronic den- 
sity and Tfi (r) are shown in the contour plots of (c) and (d) 
respectively. 



The density of states in the lUP and CP for v = 1 and 
All = 0,B = 10 T is shown in Fig. [Ill In the lUP, 
the three peaks correspond to the energies i?° (see Eq. 
([T0|) of the three levels n = 0, 1, 2 renormalized by the 
exchange interaction. The band structure of the crystal 
has three peaks, as expected for a crystal with filling fac- 
tor = 1 , but slightly displaced in energy and broadened 
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0.08 



0.016 



ql 



FIG. 9: (Color online) Effective Hartree-Fock interaction in 



due to the finite bandwidtli of each band in the crystal 
state. The lowest-energy band is fully filled while the 
other two bands are empty. From this figure, we see that 
the electron-hole continuum in the crystal state occurs in 
the energy range E G [0.6,1.2]e'^ / k£ approximately. 



B. Collective excitations, absorption spectrum and 
quantum Hall plateaus 

The absorption spectrum for the lUP has one peak at 
the frequency ujiup (q = 0) calculated in Eq. which 
corresponds to the gap in the first collective mode. The 
second dispersive mode shows up in the response func- 
tions X' fd,2) (0,2) and since are not part of the dipole 

Px,y iPx,y 

definition, it does not lead to electromagnetic absorption. 

The collective mode spectrum in the crystal phase is 
more complex. Fig. [T2l shows the dispersion relations for 
i7 = 1, i? = 10 T and A^^ = 0. Only the first low-energy 
modes that are below the electron-hole continuum are 
shown. The first, gapless, mode is the magnetophonon 
mode. It has the typical u) ^ q^l^ dispersion associated 
with the magnetophonon mode of a Wigner crystal^^. 
The other modes are gapped and correspond to more 
local deformations of the density. All modes are accom- 
panied by fluctuations of the pseudospins. 

We can get an idea of the nature of the mode 
at q = by computing the response functions 
(q = 0, q = 0,aj) with a = a;, y and = 



(0, 1) , (1, 2) , (0, 2) from the pseudospins deflned in Sec. 
IV. We find that the first gapped mode appears as a 



0.06 - 



CM 

(U 



0.04 - 



liJ 



0.02 - 



0.00 




U.01 



0.02 



0.03 



0.04 



p.ooo 



AB/(e'/(Kl) 



-0.62 
-0.64 
-0.66 



"a, -0.68 



LU 



-0.70 
-0.72 
-0.74 h 



lUP 





Fock (B) 
Fock (T) 
Total (B) 
Total (T) 



(b) 



0.01 



0.02 



0.03 



0.04 



A,/(e^/(Kl) 



FIG. 10: (Color online) Contribution of different terms to the 
total energy of the crystal state in bilayer (B) and trilayer 
(T) graphene for v = \ and i? = 10 T. (a) Hartree and bias 
energies; (b) Fock and total energies. The dashed line is the 
total energy of the incoherent uniform phase. 



pole of x^fo^.i) (0,1) and x^fi-',2) d^^, but not of x^fo^,^) (0,2) 

Pa; iPx Px iPx Px iPx 

while it is just the opposite for the second gapped mode. 
We can thus expect that the first gapped mode will be 
active in the absorption while the second gapped mode 
will not. This is confirmed by a direct calculation of the 
absorption spectrum (w) as shown in Fig. [13] for the 
crystal state at A^i = 0.09e^/K£. Fig. [T^ also shows 
the change in the absorption spectrum when is in- 
creased from All = O.OQc^/k^, in the crystal phase, to 
All = O.le^/K^ where the C2DEG has transited to the 
incoherent uniform phase. Such a change should be ob- 
servable experimentally. 

The Coulomb energy j kI = 56.2-\/B meV with B 
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FIG. 11: (Color online) Density of states in the incoherent 
uniform phase (lUP) and coherent crystal phase (CCP) for 
/7 = 1, All = and B = 10 T. 



^ 0.3 - 




lU 0.2 - 



0.5 1.0 

k/(27t/a) 



FIG. 13: (Color online) Comparison of the absorption in the 
incoherent uniform phase at All ~ O-Ic^/kI (dashed blue 
line) and in the coherent crystal phase (full black line) at 
All = 0.09eV«^- We choose i7 = 1 and B = 10 T in both 
cases. 



be a corresponding absorption at the pinning frequency 
corresponding to this gap. The pinning frequency would 
then depend on the level of disorder in the system. 

The transition from the lUP to CP should also show 
up in the Hall conductivity. As we just said, the crystal 
state would be pinned by disorder leading to insulating 
behavior of the electrons in the uppermost Landau level. 
When 1/ — 1,2, this means that the quantized Hall con- 
ductivity should have a value corresponding to the ad- 
jacent interaction driven integer quantum Hall plateau. 
For example, if the ground state stays crystalline in a 
small region around filling factor v = —5 (corresponding 
to ? = 1), there should be a dip in the Hall plateau from 
i^xij = —he^/h to Uxy — —6e^/h in this region. 



IX. CONCLUSION 



FIG. 12: Dispersion relation of the collective modes of the 
crystal phase at = 1, All ~ and B = 10 T. 



in Tesla and k = l.At B = 10 T, this gives e'^/nl = 
4.3 X 10^"^ Hz. The frequency of the first gapped mode 
is at the upper limit of the microwave spectrum. It can 
be pushed down, however, by increasing the dielectric 
constant k of the substrate. Note that the crystal state 
would likely be pinned by disorder. In this case, the 
magnetophonon mode would be gapped and there would 



We have studied the phase diagram of the two- 
dimensional electron gas in an ABC-stacked graphene tri- 
layer in Landau level iV = 0. Our analysis is restricted to 
the integer filling factors such as = — 5, — 4 and z/ = 4, 5 
where spin and layer pseudospin degrees of freedom can 
be considered as frozen. We used a three-level system 
consisting of three Landau orbitals n = 0,1,2 which are 
separated by an energy gap All (this gap is related to 
an applied potential bias between the outermost layers) 
and considered Coulomb interaction in the Hartree-Fock 
approximation. 

We calculated the energies of different uniform and 
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non-uniform phases of the C2DEG as A^^ is varied 
from negative to positive values at integer fihing fac- 
tors V — 1,2 of the three- level system. Our results show 
a phase diagram for the ABC-stacked trilayer graphene 
that is very different from its bilayer cousin. The ground 
state in the trilayer is an incoherent uniform phase at 
large value of IA^lI and a coherent crystal phase in be- 
tween. In the incoherent uniform phase, electrons fully 
occupy the lowest-energy orbitals. In the coherent crystal 
phase, there is one electron per unit cell of a triangular 
crystal lattice and this electron is in a linear combina- 
tion of the three orbital states. Around each lattice site 
of the crystal, we find an intricate orbital pseudospin tex- 
ture with a finite density of electric dipoles in the plane 
of the layers. The fact that the crystal state also occurs 
at positive value of A^^ is specific to trilayer graphene 
or, more generally, to C2DEG with chirality index J > 2. 

We studied some transport and optical properties of 
the two phases involved in the phase diagram. The crys- 
tal has a phonon mode that, when disorder is included 
in the calculation, should be gapped at q = 0. The cor- 
responding pinning mode should then be observable in 
microwave absorption experiment. The crystal has also 
higher-energy modes which are active in electromagnetic 
absorption. By contrast, only one mode should show 
up in the absorption spectrum of the incoherent uniform 
phase. It should thus be possible to localize the transi- 
tion between these two states experimentally. 

Another experimental signature of the crystal state is 
that, if disorder localizes all electrons in the uppermost 
Landau level, then a measurement of the Hall conductiv- 
ity should show a quantization at a value corresponding 
to the adjacent interaction driven integer quantum Hall 
plateau. 
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Appendix B: HARTREE-FOCK AND GRPA 
EQUATIONS 

The equation of motion for the Green's function in the 
Matsubara formalism and in the Hartree-Fock approxi- 
mation is given by 

Til. ,122 {q,iujn) (Bl) 

- "^q-q'^nf "3 (q - q') G'„3,«2 (q, i^n) 
q' 

+ Tq.q'^nf?"3 (q ~ q') '^"3 ,"2 (q, ^^n) 
q' 

— <^ni,n2^q,0) 



Appendix A: FOCK INTERACTIONS IN THE 
UNIFORM PHASES 



where 



7q,. 



(B2) 



In the uniform phases, we need to evaluate the Hartree 
and Pock interactions for q = 0. The Hartree interactions 
are zero and the only nonzero Fock interactions are (in 
units of y^n/2e'^/K£): 



^0,0,0,0 (0)-l,Xi,i,i,i(0) = -, 



Aa 



(0) 



41 
64' 



(Al) 



(A2) 



and 



"3, "4 

C^n3,ru (q) = ^«i,n4,«3,n2 (-q) (Pni,«2 (q)> ■ (B4) 



(-q)(pni,„2(q)>, (B3) 



The self-consistent Eq. (jBip can be put in a matrix 
form by defining super-indices. It must then be solved 
numerically in an iterative way in order to get the order 
parameters in the different orbital phases. 

In the GRPA, Xni,n2,n3,n4 (q^q';''') is the solution of 
the equation 
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(B5) 



ris ,ne ,n7 ,ns (Q ) X"? ,"8 ,"3 ,"4 i1 i 



-| xi?,„2,n5,«6 (q. q"; «^^») 



q" 

^^"5, "8, "7, "6 (q ) Xn7,"8^"3/'l4 (q 1 q : *^n) 7 

where r2„ is a bosonic Matsubura frequency and 
the Hartree-Fock two-particle Green's function 

(0) 

Xni,«2,n3,n4 



(q,q';«^^n) is given by 



[^hn^ - {El - El)] xl"U„„3.«4 (q' q'' (bs) 

%q,q' (Pni,n4 (q ^ q')> '5n2,«3 
-%q,q' (Pn3,n2 (q ^ q')) Sm.m 

-Y.<^"U^l (q - q") xSL,n3,n4 (q"' q'^ ^n) 

q" 
q" 

+ E ^Iq-'^i':^, (q - q") xL°!„„„3,«4 (q"' q'- ^n) 

q" 

~ Z] 7q,q"^n2,m (q ~ q ) Xni^,m,n3 ,n4 (q jqj^n)- 



Note that the response functions depend only on the or- 
der parameters {pn.m (q)) computed in the HFA. Eqs. 
(IB5IB6|) can be solved numerically by writing them in a 
matrix form defining super-indices. 



Appendix C: INFINITESIMAL GENERATORS OF 
SU(3) 

Our system has a SU(3) representation. The infinites- 
imal generators of this representation are given by the 
traceless Hermitian matrices 



where, in the basis (0, 1, 2): 



Aa 

2 



(CI) 



1 \ -i 

Ai = I 1 ,A2 = i I , (C2) 

000/ looo 



/l 

-1 I ,A4 
VO 0^ 

/ 




\ i 




(C3) 
(C4) 



and 





(C5) 



With these generators, we can define eight real fields 
Fa (r) = $1' (r) Ta^ (r) with the Fourier transform 

Fa (q) = J dre-'i '-c&t (r) (r) . (C6) 



These fields are related to the pseudospin fields by the 
relations: 



Fi (q) = (q) , F, (q) = pj^''^ (q) , (C7) 

F3 (q) = (q) , (C8) 

F4 (q) = 13 (q) , F, (q) = /3,pf 2) (q) , (C9) 

F, (q) - ^ ^ ^^^(1,2) ^ (^10) 



and 



Fs (q) = (po.o (q) + (q) - 2p2.2 (q)) , (Cll) 



where we have defined 





1 , 








- Pi.i) > 




= 2 ~ 





(C12) 
(C13) 
(C14) 
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